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INTRODUCTION 


The  accuracies  of  shipboard  navigational  tecliniques,  while  being  quite  adequate 
for  maritime  purposes,  are  not  sufficiently  accurate  for  scientific  evaluation  of 
data  recorded  do\Mi-rangc.  Li  oitlcr  to  correlate  tiiese  data  witli  data  obtained 
in  an  eailli -fixed  reference  system,  it  is  necessaiy  to  accurately  know  tlie  lo¬ 
cation  of  the  ship  relative  to  tliis  cartli  reference  system.  Witli  tliis  information,- 
tlic  impact  of  a  ballistic  missile  or  tlie  trajectory  of  a  space  probe  can  be  de¬ 
termined  in  tlie  same  reference  frame  as  tliat  used  by  a  fixed  land-based  meas¬ 
urement  device.  This  paper  describes  a  mathematical  method  for  accurately 
determining  the  location  of  tlie  DAMP  ship  in  an  earth  reference  system  using 
observed  tracking  data  on  a  tai^ct  whose  orbital  elements  are  known. 

•  Ill  essence,  tills  metliod  relies  on  three  basic  considerations: 

1.  The  teclinique  requires  tliat  tlic  ship  track  a  target  whose  orbital  elements 
are  known.  These  elements  will  be  considered  as  tlic  standard  orbital 
elements. 

2.  Assuming  any  ship's  locatioi  and  using  these  standard  orbi'al  elements, 
a  set  of  radar  parameters  (range,  azimuth,  elevation)  can  be  computed. 

3.  The  actual  ship's  location  will  then  be  that  location  at  which  the  computed 
radar  parameters  best  fit  the  observed  radar  parameters. 


PRESENTATION 


If  a  target  whose  orbital  elements  are  known  is  tracked  using  a  noise  and  bias 
free  radar,  the  location  of  the  radar  could  be  determined  purely  by  geometric 
considerations.  However,  since  tliis  is  not  possible;  the  more  soj[*isticated 
technique  described  here  is  required. 

* 

The  values  of  the  observed  range,  azimuth,  and  elevation  radar  data  (R^,  X  ,  E^ 

in  an  earth-fixed  system  are  assumed  to  be  normally  distributed  about  the  values 

that  would  be  observed  \vith  a  noise  and  bias  free  ra^r  (R  ,  A  ,  E  ).  It  is 

m  ^m  m'- 

further  assumed  that  these  distributions  are  independent  of  each  other  \vithin  the 
time  span  of  the  observations.  Based  on  these  assumptiops,  the  probability, 
density  of  observing  any  combination’  of  these  parameters  at  any  observation  i,  is: 
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Hereo- o-^,  and  o-^,  are  the  standard  deviations  in  range,  azimuth,  and  eleva¬ 
tion,  respectively. 


If  the  location  of  the  ship  is  determined  such  that  the  values  of  (A  ),  and 

(Em^  yield  a  maximum  probability  density  ot  having  observed  (R  ),  (A  ),  and 
.  .00^ 

(Eq)  over  all  observations,  then  the  requirements  of  coi;side.rati(Mi  3  stated  in 

the  Introduction  are  fulfilled. 
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This  probability  density  for  N  number  of  observations  becomes: 
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.  The  mean  radar  parameters  are  functicois  of  the  ship  location  and,  therefore, 

may  be  individually  expanded  about  some  arbitrary  location  c.  The  linear  terms 

•  • 
of  a*Tayl9r  series  are  used:  ^ 


8R.  SR, 
^A.  •  3A. 


in  which  \  and  L  represent  latitude  and  longitude,  respectively.  The  "c"  sub¬ 
script  indicates  evaluation  at  point  c  and  the  delta,  expression  represents  the  dif¬ 
ference  in  position  between  c  and  the  actual  ship  location. 

•  •  . 

Substituting  these  expressions  for  the  mean  radar  parameters  into  Equation  (1) 

yields: 


•  ,  • 

Since  it  is  desired  to  maximize  P.,(R  ,  A  ,  E  )  above,  then  the  exponent 

•  n  o  *  0  0  • 

should  be  miniiplzed  \vilh  respect  to  and  AL.  Therefore, 
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It  is  now  necessary  to  obtain  expresslwis  for  range,  azimuth,  and  elevation  in 
terms  of  the  latitude  and  longitude  in  order  that  the  partial  derlvitlves  appearing 
in  Relations  (2)^d  (3)  may  be  evaluated.  These  expressions.are  glOen,by^ 

*  •  *  5 

R  =  r(Rp  +  h)^  +  RE^-2RE(RE+h)  {cosAscosAt  co3(L-p-Ls) +3inAx  sin  As|]  ^ 
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where  is  the  radius  of  the  earth,  h  is  tlie  altitude  of  the  target,  K  g  and  JLjg 
are  the  latitude  and  longitude  of  the  ship  at  locatioi  c,  and  are  the  laU 
tude  and  longi^de  of  the  target. 

»  • 

The  desired  partial  derivatives.are  now,. 
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in  which 


Tliese  expressions  for  the  partial  derivatives  may  now  he  substituted  into  Rela¬ 
tions  (2)  and  (3)  and  a  soluticm  for  AL  and  obtained  based  on  all  observations. 
However,  these  values  represent  a  correction  to  the  assumed  location  aixl  may 
still  contain  error.  This  error  results  from  using  only  the  linear  terms  of 
the  Taylor  expansion. 

Successive  iterations  using  these  corrected  values  as  the  new  locatimis  yields 
a  sequence  of  values  for  AL  and  AX  \sbich  ccmvergo  toward  zero. 

•  • 

A  possibility  exists  that  a  noi-ccaivergent  sequence  may  result.  However,  in 

these  investigations,  which  assumed  rather  large  errors  in  initial  locattcm,  this 
did  not  occur.  Related  investigations  using  these  techniques  have  shown  that 
under  tliese  conditions  the  simple. expedient  of  selecting  a  different  locaticm  pro- 
duced  convergent  sequences.  *  *  • 

a 

•  *  *  • 

A  series  of  simulated  tests  were  performed  using  orbital  elements  from  an 
actual  ballistic  missile  trajectory  to  determine  the  feasibility  of  thi^  method. 

The  simulations  involved  the*following  basic  procedure. 

1.  A  location  was  chosen  from  which  a  noise-free  set  of  radar  parameters 
•  * 

were  computed.  .  •  .  . 

2.  These  parameters  were  then  contaminated  with  gaussian  noise  of  known  • 
•  •  • 

standard  deviation  to  simulated  observed  radar  parameters.. 

•  ^  •  • 

•  3.  Using  these  parameters  and  assuming  an  arbitrary  ship's  locatiwi,  in 

conjunction  ivith  the  known  orbital,  elements,  the  method  was  applied  to 

determine  the  simulated  actual  location.  •  * .  »i 

•  • 

t 

It  should  be  qoted  that  one  of  the  inputs. is  the  standard  deviation  of  the  noise, 
which  in  an  actual  mission  may  not  be  known.  In  order  to  be  realistic,  the  input 
values  of  standard  deviations  were  selected  to  be  different  from  those  used  in 

establishing  the  simulated  observed  radar  parameters. 

•  • 

•  The  values  of  the  input  standard  deviations  assigned  to  each  of  the  radar  param- 
eters  (range,  'a*zimufli,  and  elevation)  are  measures  of  their  confidence  level: 
i.  e. ,  if  one  of  the  parameters  is  heavily  contaminated  with  noise,  its  confidence 
level  is  low.  '  . 


RESULTS  AND  CONCLUSIONS 

A  table  listing  tlie  conditions  and  results  of  tlie  simulation.progi’ain  follows. 
True  ship's  location  was  selected  on  the  basio  of  a  "worst  case"  situation 

where  tlie  ship  lies  very  close  to  Uie  plane  of  Uie  trajectory. 

*  • 

Di  tlie  first  simulation  listed,  the  uiput  values  of  Uie  standard  deviation  were  Uie 
same  as  Uie  values  used  to  contaminale  Uic  radar  parameters.  Even  with  the 
raUier  large  errors  in  assumed  ship's  location,  Uie  technique  indicated  true 
location  to  slightly  more  Uian  six  feet.*.* 

*  • 

The  ne.xt  four  simulations  tested  the  method'for  varying  degrees  r,l  assumed  noise 
contaiiiination.  Note  Uiat  regardless  of  Uie  aniomit  of  error  in  Uie  assumed 
noise,  or  wjieUier  it  exceeds^or  falls  below'  Uie  true  standard  deviation,  the  lo¬ 
cation  is  determined  to.Uie  same  degree  cf  accuraey,  30  feet  in  latitude  and  241 

feet  m  longitude.  It  must  be  kept  in  mind.Uiat  Uie  tecluiique,  in  addition  to  cor¬ 
recting  for  ship's  location,  yields  corrected  information  on  standard  deviation 
Uiat  may  be  re-inserted  to  pi’ovide  a  ship's  location  wiUi  an  accuraey  commen¬ 
surate  with  that  detei-mined- in  the  first  simulatioi..  . 

*  •  • 

•  • 

•  •  •  • 

The  last  simulation  not  ondy  introduced  noise,  but  also  an  amount  of  bias  equal 

,  •  • 

to  that  of  the  noise.  Undei*  these  condittens,  -loeation  in  latitude  was  determined 

to  within  1471  feet  and  in  longitude  to.wiUiin  1050  feet.  These  results  indicate 

•  •  • 

that  bias  errors  in  the  radar  parameters  .effect  Uie  accuracy  of  determining  the 
position  of  Uie  ship.  •  . .  • 

In  any  radar  system  there  exists  many  sourees  whieh  introduee  bias  errors.* 
Knowing  these  sources,^  it  is  possible  to  extend  this  differential  eorrection  teeh- 
nique  to  eompute  and  eompensate  for  these  errors.  For  the  DAMP  tracking 
radar  the  major  bias  errors  are  introdueed  in  the  determination  of  true  north 
and  the  loeal  vertieal.  However,  ip  the  above  simulation  these  errors  were  not 
eonsidered  sinee  this  investigation  was  undertaken  primarily  to  determine  the 
feasibility  of  this  method.  The  teehnique  is  currently  being  extended  to  inelude 
bias  errors  in  true  north'and  local  vertical. 
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TABULATED  REI 
SHIP'S  LOCATION  DIF! 

COm^ECTIO] 


TIME  OF  INITIAL  RADAR 

OBSERVATION  (SECONDS  AFTER  MIDNIGHT):  53826.0  *TRIJ 

TIME  BETWEEN  OBSERVATIONS:  1.0  SECOND 
TIME  OF  FINAL  RADAR 

OBSERVATION  (SECONDS  AFTER  MIDNIGHT):  54026.  OO 


SIMULATION 

NUAIBER 

RADAR  BIAS 
ERRORS 
ASSIGNED 

TRUE  STANDARD 
.DEVIATION  . 

ASSUMED 

STANDARD  DEVIATION 

ASSUMED  SHI 
POSITION 

(Degrees) 

(Yards) 

(Degrees) 

(Yards) 

(Degrees) 

(Yards) 

(Degrees) 

A 

E 

R 

A 

E 

R 

A 

E 

R 

Latitude 

Long 

1 

• 

0 

0 

0 

.0.25, 

0.25 

15.0 

0.25 

0. 25 

15.0 

8.713972S 

17. 

• 

2 

0 

0 

.0 

0.25 

0.25 

15.0 

0.5  , 

0.5 

15.0 

8.713972S 

17. 

3 

0  . 

0 

0 

0. 25 

0. 25 

15.0 

2.5 

2.5 

15.0 

8.713972S 

17. 

4 

0 

0 

p 

0.25 

• 

• 

0.25 

15.0 

0.125 

0.125* 

15.0 

• 

8.713972S, 

•17. 

5  ■  . 

0.25 

0. 25 

15.0 

25  X  10® 

25  X  10® 

15.0 

8.713972S 

17. 

6 

• 

0.25 

0.25 

15.0 

0. 25 

0.  25 

15.0 

0.25 

0.2^ 

15.0 

8.713972S 

17. 

TABULATED  RESULTS 
SHIP’S  LOCATION  DIFFERENTIAL 
CORRECTION 

J26.0  TRUE  SHIP'S  POSITION:  7.713972*  SOUTH  , 

'  •  .  19.00*  WEST 


}26. 00 


ASSUMED 

\NDARD  DEVIATION 
Degrees)  (Yards) 

ASSUMED  SHIP'S 
POSITION 
(Degrees) 

MINIMUM 
NUMBER  OF 
ITERATIONS 
REQUIRED 

FIRST  ITERATION 
POSITION 
(Degrees)  * 

LAST  ITERATION 
POSITION 
(Degrees) 

ERROR  IN  LAST 
ITERATION  (FEET) 
FROM  TRUE  SHIP'S 
POSITION 

,  • 

E 

R 

Latitude 

Longitude 

Latitude 

Longitude 

Latitude 

Longitude 

Latitude 

Longitude 

0.25 

15.0 

8.713972S 

17.  OW 

4 

8.  248942S 

19.41574W 

7.713993S 

19. 00002W 

6.6 

6.3 

• 

0.5 

15.0 

8.713972S 

17.  OW 

4 

• 

• 

8. 249214S 

19.41537W 

7.713876S 

18.99923W 

30.4 

241 

2.5 

15.0' 

• 

8.713972S 

.  17.  OW 

4 

8. 249509S 

19.4151 

>4W 

7.713877S 

18.99924W 

30 

241 

0.125 

15.0 

8.713972S 

17.  OW 

4 

8*244584S 

19.4in9W 

7.713876S 

18.99924W 

30.4 

241 

0« 

25  X  10® 

15.0 

8.713972S 

17.  OW 

4 

« 

8.  249529S 

19.41566W 

7^7138748 

18.99924W 

29.7 

241 

0.25 

15.0 

8.713972S 

■ 

17.  OW 

4 

8.240465S* 

19.409.19W 

7.709324S 

18.99668W 

1471 

1050 

TRAJECTORY  ORBITAL  ELEMENTS  . 

#  ' 

Semi-major  axis:  0.82363  Argument  of  Perigee :  335.4405 

Time  of  Perigee:  51464. 27  Eccentricity:  0. 43272 

Right  Ascension:  25.27468  Inclination:  31.80759 
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